Gapped spin liquid phase in the J\ — Ji Heisenberg model by a Bosonic resonating 

valence-bond ansatz 
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We study the ground-state phase diagram of the spin-1/2 J\— J2 Heisenberg model on the square 
lattice with an accurate Bosonic resonating valence-bond (RVB) wave function. In contrast to the 
RVB ansatz based on Schwinger Fermions, the representation based on Schwinger Bosons, supple- 
mented by a variational Monte Carlo technique enforcing the exact projection onto the physical 
subspace, is able to describe a fully gapped spin liquid in the strongly frustrated regime. In partic- 
ular, a fully symmetric Z2 spin liquid is stable between two antiferromagnetic phases; a continuous 
transition at J2 = 0.4Ji, when the Marshall sign rule begins to be essentially violated, and a first- 
order transition around J2 = O.6J1 are present. Most importantly, the triplet gap is found to 
have a non-monotonic behavior, reaching a maximum around J2 = 0.51Ji, when the lowest spinon 
excitation moves from the F to the M point, i.e., k = (n,0). 
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I. INTRODUCTION 

The search for quantum spin liquids in frustrated 
quantum antiferromagnets has a long history* 1 - In recent 
years, thanks to the advance of numerical techniques, 
several candidates for spin liquids have emerged in two- 
dimensional (2D) systems. These include the Hubbard 
model on the honeycomb lattice^ the spin-1/2 Heisen- 
berg model on Kagome lattice^ and more recently the 
spin-1/2 Ji— J 2 Heisenberg model on the square lat- 
ticed In all these cases, a small but finite spin gap has 
been found and, according to generalizations of the Lieb- 
Schultz-Mattis theorem for higher dimensionalities ^ a 
topological degeneracy is expected. In spite of these 
results, descriptions based upon Fermionic resonating 
valence-bond (RVB) theory predict more often the ex- 
istence of gapless spin-liquid states. For example, for 
the J1 — J2 model on the square lattice^ the Heisenberg 
model on the triangular— or the Kagome lattices^ and, 
more recently, also for the unfrustrated honeycomb lat- 
tice^ the Fermionic RVB theory always predicts a gapless 
spin liquid phase with a Dirac-type spinon dispersion as 
the best variational state. 

The J\— J 2 model represents the simplest model to 
study the effect of frustration in a (low-dimensional) 
magnetic system; for this reason it has been investigated 
by many different approaches in the last 20 yearsj At 
the classical level, the system is magnetically ordered for 
J 2 < 0.5Ji with the standard antiferromagnetic pattern 
at q = (7T,7r). For J 2 > 0.5Ji, the ordering wave vector 
is moved to q = (ir, 0) or (0, 7r); these two ordered phases 
are separated by a first-order transition. Within the lin- 
ear spin-wave approach, which goes beyond the classical 
theory, quantum fluctuations destroy the magnetic or- 
der in the intermediate region of 0.4Ji < J 2 < O.6J1, 
hence leading to a magnetically disordered stated 7 How- 
ever, the nature of this disordered phase is still elu- 



sive and several proposals have been raised. These in- 
clude valence-bond solids with broken spatial symme- 
trie a 11 ' 14 ' 18 ' 19 or gapless spin-liquid states^ The latter 
proposal is especially attractive, since it provides a sim- 
ple and very accurate Fermionic RVB wave function for 
0.4Ji < J 2 < 0.55Ji. This state has a Dirac-type spinon 
dispersion and Z 2 gauge structure and becomes stable for 
J 2 > 0.4 Ji. 

More recently, density-matrix renormalization group 
(DMRG) calculations provided some evidence for a fully 
gapped spin liquid in the intermediate region of 0.4 J] < 



J~2 ^5 0.62Ji. 4 Within this numerical approach, the spin 
gap increases linearly from J 2 ~ 0.4Ji, reaches a max- 
imum around J 2 ~ 0.59Ji, and then rapidly decreases. 
For J 2 > 0.62Ji, a collinear magnetic order takes place. 
The spin-liquid phase determined by these DMRG cal- 
culations is thus inconsistent with the Fermionic RVB 
theory, due to presence of a finite spin gap. 

In this paper, we investigate the spin-liquid phase of 
the J\ — J 2 model with a Bosonic RVB wave function. 20 
This is motivated by the following reasons. First, while 
the Fermionic RVB state is found to be unable to open 
a spin gap for this system, a Bosonic spin-liquid state 
is by definition gapped, because otherwise the (Bosonic) 
spinon would condense and the system would develop 
magnetic order. Second, since the spin-liquid phase is 
found to exist in a quite small region between two mag- 
netically ordered phases (for which a Bosonic description 
is quite accurate), it is natural to expect that the inter- 
mediate spin-liquid phase inherits some Bosonic charac- 
teristic. 

The Bosonic RVB state has been adopted in many pre- 
vious studie a 21 i 22 and is found to describe quite well both 
the magnetic ordered state and the disordered state for 
unfrustrated systems i 23 ' 24 For frustrated magnetic sys- 
tems, the use of the Bosonic RVB wave function is very 
limited, since the loop gas algorithm for the Bosonic RVB 
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state encounters serious sign problems; moreover, the 
computation of the wave function amplitude in the or- 
thogonal Ising basis involves permanents of matrices,— 
implying a computational cost that grows exponentially 
with the size of the system. Only very recently, this ap- 
proach has been implemented on small clusters for the 
Kagome latticed 

Here, the Bosonic RVB state is obtained after project- 
ing the ground state of the mean-field Schwinger Bo- 
son Hamiltonian^T. into the physical subspace with one 
spin per site. After this projection, the wave func- 
tion turns out to be equivalent to the standard Liang- 
Doucot- Anderson RVB ansatz, 20 defined only in terms of 
a bosonic pairing function (that connects opposite sub- 
lattices) . To enforce the physical symmetry of the model 
in the RVB state, we have made a full symmetry clas- 
sification of the Schwinger Boson mean-field ansatz on 
the square lattice with the projective symmetry group 
(PSG) technique.— ~— Then, we have performed varia- 
tional Monte Carlo simulations in order to optimize such 
a Bosonic RVB state, by using both the permanent Monte 
Carlo algorithm and the loop gas algorithm. 

We find that the Bosonic RVB wave function gives a 
rather good variational description of the system. In 
addition, we find that the phase diagram predicted by 
the DMRG calculations can be well reproduced. More 
specifically, the system is found to enter a fully gapped 
spin liquid state around J2 = 0.4Ji through a continu- 
ous transition, when the Marshall sign rule in the ground 
state begins to be essentially violated. A level crossing 
of the spinon excitation is observed around J2 = 0.51Ji, 
when the gap minimum of the spinon excitation branch 
is moved from the V to the M (i.e., k = (tt, 0)) point and 
a kink appears in the spin gap as a function of J2. 

Finally, by PSG symmetry considerations, it can be 
shown that the spin gap is always finite at the M point 
in the spin-liquid region (while it can vanish at (tt, tt), at 
the transition to the antiferromagnetic phase for small 
J2). This fact implies that the magnetic structure factor 
is always finite at the M point, ruling out a continuous 
transition to the collinear phase at large J 2 - 

The paper is organized as follows: in Sec. HH we de- 
scribe the model and the method; in Sec. IIII1 we present 
our numerical results; finally, in Sec. IIV| we draw our 
conclusions. 

II. THE MODEL AND METHODS 

In this paper, we consider the following model: 

H= J 1 ^S i -S i + J 2 ]T S.-S,, (1) 

<(i,i» 

where and indicate nearest-neighbor and 

next-nearest-neighbor sites on the square lattice, respec- 
tively; Si denotes the spin operator at site i. 

In the Schwinger Boson representation^ the spin op- 
erator is written as S = \ J2 a 8 ba&a,pbp, where b a is a 



Boson operator, a is the Pauli matrix. Bosons should sat- 
isfy the no double occupancy constraint J2 a b^ a b a = 1, in 
order to be a faithful representation of the spin-1/2 oper- 
ator. Within this representation, the Heisenberg super- 
exchange coupling can be written as (apart from addi- 
tive constants) Sj • Sj = —^A\jAi t j — \B\jBi.j, where 
= hfbji - b ti b n and Bij = b\^b n + b^b^ZL 

In the mean-field treatment, we replace Ajj and Bij 
with their mean-field expectation value Aij and Bij, so 
to have: 

Hmf — - -Ys^ijAlj + h.c.) 

i,j 
i,3 

+ a e(e & l^-iY (2) 

where Ajj = JijAij, Fij = JijBij, and the chemical 
potential A is introduced to fulfill, on average, the single- 
occupancy constraint. The mean-field ground state has 
the general form of 

|MF) oc exp j^a^i^X^ - &j^t f ) J |0). (3) 

Then, a suitable RVB wave function in the physical 
Hilbert space with one Boson per site may be obtained 
by projecting the mean-field state, namely 

|RVB) =P G |MF), (4) 

where Pq is a Gutzwiller projector that enforces the con- 
straint of one Boson per site. The equivalence of the RVB 
state with the standard Liang-Doucot- Anderson stated 
is clear after projection onto the physical subspace. 

The form of the RVB amplitude a(Ri,Rj) is deter- 
mined by the parameters Ajj, Fi t j and A. At the mean- 
field level, Ajj and F^j are non-zero only on those bonds 
with Ji j =/= 0. However, from the variational point of 
view, we can take {Ajj, Fij, A} as a set of free parame- 
ters to construct the RVB state. In such a case, we can 
also introduce Ajj and Fi j on longer bonds, for which 
Ji,j = 0. 

In order to describe a spin liquid state with the 
full symmetry of the model, the mean-filed parameters 
{Aij,Fij,X} must satisfy certain symmetry conditions. 
Since there exists a U(l) gauge degree of freedom in the 
Schwinger Boson representation of the spin operator (i.e., 
bi, a &j, a e'" leaves Si unchanged), the symmetry re- 
quirement on the mean-field Hamiltonian is actually the 
U(l) gauge projective extension of the physical symmetry 
of the model. Such symmetry conditions on the mean- 
field ansatz can be readily worked out by the so called 
PSG technique developed by Wen 28 for the Fermionic 
representation. The Bosonic version of the PSG is the 
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U(l) subset of the Fermionic PSG. 29,30 Here, we will just 
point out some basic structures that are relevant to our 
study. 

In the Schwinger Boson formalism, the mean-field pa- 
rameters Aij and Fij describe antiferromagnetic and 
ferromagnetic local correlations, respectively (see Ap- 
pendix [A] for the possible phases implied by this ansatz). 
Here, we assume a non-zero Ajj between nearest- 
neighbor sites. Then, we find that a non-zero Ajj be- 
tween next-nearest-neighbor sites is compatible only with 
the so-called type B translational property of the mean- 
field Hamiltonian ) 28 ' 29 which implies a unit cell with two 
sites. We find that such state is much higher in energy 
than any state in the so-called type A class, characterized 
by a manifestly translational invariant mean-field ansatz. 
Therefore, in the following we restrict our analysis only 
to translationally invariant states. Within the type A 
states, we find the following general rules for the mean- 
field ansatz for a symmetric spin liquid state. First, for 
sites belonging to different sub-lattices, only a real Ajj 
is allowed. Second, for sites in the same sub-lattice, only 
a real Fij is allowed. Considering the site i as belonging 
to A sub-lattice^ the allowed mean-field parameters up 
to the fourth-neighbor are given by: 
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where <5 M (with fi = 1, . . . , 4) denotes the vectors connect- 
ing the site i to its neighbors, up to the fourth distance. 
Here, {\, F, A. F 2x , A 2xy } are a set of real parameters. 
An illustration of the ansatz used in this study is shown 
in Fig. [U For the sites i belonging to B sub-lattice, the 
sign of A and A 2xy should be reversed (since Ajj is odd 
by interchanging i and j). 

At the mean-field level, both F 2x and A 2xy are zero, 
and the Hamiltonian is given by 
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FIG. 1: (Color on-line) An illustration of the mean-field pa- 
rameters starting from a site in the sub-lattice A. Gray and 
dark dots denote sites in sub-lattices A and B. Here, only 
bonds up to the third neighbors are reported, since longer- 
range parameters are found to be negligibly small after opti- 
mization. The pairing term is always directed from sub-lattice 
A to sub-lattice B 




in which MBZ indicates the reduced (magnetic) Bril- 
louin zone, tpl = (&\ kt , &^ kf , 6^-14, ^B-kj.), e k = A + 
2Fg(k), and A k = 2A7(k). Here g{k) = cos^) cos(k y ), 
7(k) = (cos^) + cos(k y ))/2. The mean-field spectrum 
is given by £" k = -\/ e k ~ A k and the minimal spinon gap 
is given by 



\/(A + 2^) 2 
X-2F, 



{2Af 



2XF < A 2 ; 
2AF > A 2 . 



For the first case, the gap minimum is located at the T 
point, while for the second case the gap minimum is at 
the M point. 



FIG. 2: (Color on-line) The optimized values for various pa- 
rameters: 2F/X and 2A/A (a). The thin lines denote the 
solution of the mean-field self-consistent equations. 2F2 X /\ 
and the (normalized) spinon gap at the F point (b). 



Finally, the RVB amplitudes derived from the mean- 
field ground state are given by 



a(i?i — Rj 



1 E 



N ^ e k + £ k 

keMBZ 



(10) 



where N is the number of sites, i € A and j G B. The 
RVB amplitudes between sites in the same sub-lattice are 
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identically zero. We would like to mention that, within 
the standard formulation based upon Monte Carlo sam- 
pling ,22^22 only positive pairing functions a(Ri-Rj) have 
been considered so far. In our formulation this restric- 
tion applies only for standard antiferromagnetic phases, 
while negative amplitudes are found in the much more 
interesting spin-liquid phase. 



III. RESULTS 

The mean-field Hamiltonian ([9]) has been studied by 
Mila and collaborators, .22, showing that no spin-liquid 
phases are stabilized and a direct transition between two 
ordered phases is present, with a phase diagram that is 
very similar to the classical limit. 

In order to go beyond this approximation, we now 
move to the projected RVB state of Eq. to assess 
the possibility that quantum fluctuations may induce a 
finite spin gap and, therefore, a stable spin liquid. We 
thus determine the parameters in the Bosonic RVB state 
by optimizing the energy of the original J1—J2 model, 
rather than solving the self-consistent equations. Then, 
the spinon gap can be estimated by inserting back the 
optimized parameters into the mean-field dispersion re- 
lation E-k. Note that the RVB wave function does not 
depend on the overall energy scale of the system. As a 
result, the spinon gap can be determined only up to a 
normalization constant. Here we will use the chemical 
potential A as the unit of energy. To have an estimate 
of the absolute scale of the spinon gap, we determine the 
pairing potential A from the equation 



A = Ji{A 



i,i-\-x / 
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N 



E 

keMBZ 



A k7 (k) 

Ey 



(11) 



by inserting on the right-hand side the optimized values 
of A/A and F/X, which are (A/A) opt and (F/X) opt . Then 
A can be determined by requiring that A/ A = (A/A) op t. 

The computation of the Bosonic RVB wave function 
is very expensive in the Ising basis, since it requires the 
calculation of permanents, for which no polynomial algo- 
rithm exists. 25 However, on small clusters the calculation 
is still affordable. In this work, we have used a 6 x 6 clus- 
ter to perform the optimization of the parameters in the 
RVB wave function^ It is important to note the key 
difference between the mean-field theory and the pro- 
jected RVB wave function. In the mean-field theory, the 
chemical potential A is determined by the self-consistent 
equation for the total Boson number. When the spinon 
gap approaches zero, the number of Boson will diverge. 
Thus, on any finite lattice, the spinon gap can never be 
zero and a finite-size gap must exist (see Appendix [B] for 
the details on the spinon gap in the mean-field approach) . 
On the contrary, after projection, the constraint of one 
Boson per site is satisfied exactly and such a divergence 
will not appear. Therefore, the RVB wave function is 
well behaved even when the spinon gap is zero. This fact 



implies that a vanishing spinon gap can be realized ex- 
actly after optimization of the corresponding projected 
RVB wave function on a relatively small cluster. 

From our numerical optimization, we find that a spin 
gap can not be opened if we keep A/ A and F/X only. 
Moreover, by a direct optimization of the pairing ampli- 
tudes a(Ri — Rj), a good accuracy can be achieved only 
by including a third-neighbor parameter i^^/A, while 
the fourth-neighbor parameter A2 Xy /X is found to al- 
ways negligibly small. Therefore, in the following, we 
optimize the wave function with A/A, F/X and Fi x /X as 
variational parameters. In particular, we find that the 
inclusion of i^x/A is crucial for the opening of the spin 
gap. The optimized value of the parameters in the RVB 
wave function are shown in Fig. [2l 

The spinon gap at the T and the M points is shown 
in Fig. [3](a). Around J2 = 0.51Ji, a level crossing in the 
spinon excitation occurs and the gap minimum changes 
from r to M. By further increasing J2, the spinon gap 
at M decreases and eventually approaches zero around 
J2 = O.6J1. At this point the system becomes unsta- 
ble with respect to magnetic ordering at q — (tt,0). It 
should be noted that, although the spinon gap at the M 
point approaches zero continuously for J2 = O.6J1, our 
state cannot be continuously connected to the collinear 
ordered state, and a first-order transition must exist be- 
tween the fully symmetric spin liquid and the collinear 
ordered magnetic phase. This is clearly seen in the static 
spin structure factor: 



S(q) 



A k A, 



1 \ "> / £k€q-k 

2 ^ \ E k E q _ k 



q— k 



1 



(12) 



Since Aj^^o) = by symmetry (see Appendix [Cj , the 
singularity in the coherence factor for -Ek=(7r,o) — > is 
removed and the spin structure factor at q = (tt, 0) is 
always finite. Thus, the state cannot be connected to 
the collinear ordered phase, in which S(tt, 0) diverges. 
Therefore, we conclude that a first-order transition must 
exist between the spin liquid and the collinear ordered 
phase. 

Given the results for the spinon spectrum of Fig. [21a,), 
it is possible to make some prediction on the behavior of 
the triplet gap as a function of J2. Indeed, to construct 
a triplet excitation at q = (tt, tt), we can use two spinons 
both from the T point and the M points On the con- 
trary, for a triplet excitation with momentum q = (tt, 0), 
we should use one spinon from the T point and another 
spinon from the M point. Therefore, the lowest triplet 
excitation is always realized at q = (tt, tt) and the energy 
of triplet excitation at q = (tt, 0) is always finite, see 
Fig- 121b). This is consistent with the result of the static 
spin structure factor mentioned above and points to the 
fact that our spin-liquid state cannot be continuously 
connected to the collinear ordered phase. We note that 
the peculiar behavior of the triplet excitations found in 
this work represents an astonishing consequence of frac- 
tionalized spinon excitations in the spin-liquid phase. 
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FIG. 3: (Color on-line) The normalized spinon gap at T and 
M points in the spin liquid regime (a). Normalized triplet 
gap at q = (it, tt) and q = (tt, 0) (b). 
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FIG. 4: (Color on-line) Accuracy of the ground-state energy 
calculated from the best Fermionic of Ref.|6|and Bosonic RVB 
variational wave functions on a 6 x 6 lattice. 



We would like to mention that our results for the spin 
gap are quite similar to the DMRG ones<^ Indeed, within 
both approaches, the spin gap is found to open around 
Ji = 0.4 Ji and close around Ji = 0.62Ji. In addition, 
a sharp maximum is present, though its position in the 
Bosonic RVB approach is found to correspond to a lower 
value of Ji with respect to the DMRG study. Moreover, 
taking the value of A estimated from Eq.ll, which is A ~ 
1.02Ji at J2 = 0.5 Ji, we have that the maximal spin gap 
is quite consistent with the DMRG prediction. 

In this work, the sharp maximum in the spin gap is in- 
terpreted as the result of a level crossing in the minimum 



FIG. 5: (Color on-line) The normalized RVB amplitudes 
a(Ri,Rj) on different distances on a 6 x 6 lattice as func- 
tions of J2/J1. The nearest-neighbor (1,0) amplitude has 
been taken equal to one (a). The average Marshall sign of 
Eq. (|13[) calculated from the Bosonic and Fermionic RVB wave 
functions and the exact ground state on a 6 x 6 lattice as func- 
tions of J2/J1 (b). 



of the spinon spectrum (from the F to the M point). In 
such a picture the lowest triplet excitation within the 
symmetric spin liquid phase is always at q = (7r,7r). 
However, other possibilities for this structure may exist, 
among which a spin nematic liquid phase, which breaks 
the reflection symmetry x — > y but with all other phys- 
ical symmetries intact, is especially interesting. 35 Since 
the DMRG calculations have been done on rectangular 
clusters, the nematic liquid phase can be connected to 
the symmetric state continuously on finite lattices. 

To further check the accuracy of the Bosonic RVB wave 
function, we computed the relative error in the ground- 
state energy, namely AE = \E — E val \/\E Q \, where E 
is the exact ground-state energy and E val is the vari- 
ational energy of the RVB state. In Fig. |H we report 
the accuracy of the Bosonic RVB wave function on the 
6x6 cluster, in comparison with the best Fermionic RVB 
wave function.— For small J2, the Bosonic RVB wave 
function is much more accurate than the Fermionic RVB 
wave function, which cannot describe magnetically or- 
dered states. In this region, our results for the Bosonic 
wave function agree with previous calculations reported 
in Ref. 36|, obtained with a different algorithm^! or a dif- 
ferent parametrizationi22 For Ji > 0.45Ji, the Fermionic 
wave function becomes more accurate. However, the er- 
ror in both wave functions are similar and both increase 
with the same trend by increasing J 2 up to J 2 = O.6J1. 

As pointed out in Ref. [6j, the sign structure of the 
ground state is crucial for the origin of the spin liquid 
phase. For J2 = 0, the ground-state wave function satis- 
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fies the Marshall sign rulei^i However, the Marshall sign 
rule is essentially violated only for J 2 > 0.4 J-y and, in the 
Fermionic RVB approach, a Z 2 spin liquid phase emerges 
just at the same points A similar scenario also appear 
in the Bosonic representation. In this case, when the 
RVB amplitudes from sub-lattice A to sub-lattice B are 
positive, then the wave function satisfies the Marshall 
sign rule, otherwise (if some amplitudes are negative) 
the Marshall sign rule is violated. In Fig. [SJ we plot all 
the independent RVB amplitudes a(i? i ,i? J ) on a 6 x 6 
lattice of the optimized wave function (with the ampli- 
tude between the nearest-neighbor sites equal to one). 
For J 2 < 0.4 Ji, all amplitudes are positive and thus the 
wave function has the Marshall sign. For J 2 > 0.4Ji, the 
amplitude on bond (1, 2) becomes negative and the Mar- 
shall sign rule is violated. It is just at this point that the 
spin gap opens. Thus, the origin of the spin gap and the 
existence of the spin liquid phase can be understood as a 
result of violation of the Marshall sign rule. Such an un- 
derstanding is consistent with several previous studies^ 
in which the topological degeneracy, which is a hallmark 
of gapped spin liquid, is argued to be absent in system 
satisfying the Marshall sign rule. 

Finally, we report in Fig. [5] the average Marshall signs 
in the Bosonic and Fermionic RVB wave functions: 

(S) =Y / \(x\RVB)\ 2 sign{(x\RVB)(-l) N ^}, (13) 

X 

where \RVB) denotes the RVB variational state (either 
Bosonic or Fermionic) and the sum is over the orthogonal 
Ising basis \x); for comparison, we also report the results 
for the exact ground state, where \RVB) is replaced by 
|^o)- The Fermionic RVB wave function is better in the 
sense of sign structure and this is consistent with the 
fact that the Fermionic wave function has a lower en- 
ergy for large J 2 . However, it is clearly seen that both 
the Bosonic and the Fermionic RVB wave function un- 
derestimate seriously the frustration of the sign in the 
spin-liquid regime. 




Fix 



FIG. 6: (Color on-line) Phases described by the Bosonic RVB 
wave function |RVB) in the parameter space (F/X, A/ X). 



Marshall sign rule. This fact provides strong support for 
previous arguments for the absence of topological order 
in systems satisfying the Marshall sign rule. Despite that 
these outcomes are in good agreement with recent DMRG 
calculations of Ref. [H, the gapless Dirac-type Fermionic 
RVB ansatz remains slightly more accurate at the varia- 
tional level in the highly-frustrated regime. 
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IV. CONCLUSIONS 

In conclusion, we find the Bosonic RVB wave function 
generates a ground-state phase diagram of the J\ — J% 
model on the square lattice that is qualitatively consis- 
tent with DMRG results. A gapped spin-liquid phase 
is found for 0.4J! < J 2 < O.6J1. The spin-liquid 
phase is connected to the staggered magnetic ordered 
state through a continuous transition but cannot be con- 
nected continuously to the collinear magnetic ordered 
state and a first-order transition between the two must 
exist. The spin gap is found to have a maximum around 
J 2 = 0.51Ji, as a result of the level crossing between the 
spinon at T and M points. This fact implies that the low- 
est triplet excitation is found to be always at q = (jr, w) 
in the spin-liquid phase. We also found that the spin gap 
opens at the same point where the system violates the 



Appendix A: The various phases described by the 
wave function studied in this work 

The various phases described by the wave function 
studied in this work are shown in Fig. [6] Here, we report 
the various properties as a function of two parameters, 
namely F/X and A/A. The case with non-zero -F^/A is 
qualitatively similar. The condensation lines denote the 
magnetically ordered states with staggered or collinear 
patterns. The three regions, A, B and C, denote spin- 
liquid phases. In regions A and B, the spinon gap min- 
imum is realized at the T point, while in the region C 
the gap minimum is moved to the M point. In region B, 
the Marshall sign rule is violated while it is satisfied in 
region A. The region with slanted lines is physically non 
accessible. 



7 




FIG. 7: (Color on-line) The spinon gap predicted by mean- 
field theory both at F and M points. For the F point, we 
show the scaling of the mean-field gap with the linear size 
of the system. The mean-field gap at the M point in the 
thermodynamic limit is also reported (dashed line). 

Appendix B: The mean-field finite size gap 

In the mean-field theory, the spinon is always gapped 
when the system is defined on a finite lattice. We re- 
port in Fig. [7] the spinon gap obtained by solving the 
mean-field self-consistent equations. Here, we would like 
to emphasize that the origin of a spinon gap obtained 
on finite lattices with the projected Bosonic RVB wave 
function is totally different from that obtained within the 
mean-field approximation. Indeed, after projection, the 
number of spinons is fixed (each site is occupied by one 
and only one spinon) and the RVB wave function is al- 
ways well defined. 

In fact, we find that the spinon gap is exactly zero for 
Ji < QAJ\ from our optimization on the 6x6 lattice, 
see Fig. [3] Instead, the finite size gap in the mean-field 
theory is much larger and smoother than that obtained 
with the projected Bosonic RVB wave function. In addi- 
tion, we note that the mean-field theory always predicts 
a very large gap at the M point in the thermodynamic 
limit. 



Appendix C: The proof of A k= (^ m — 

For the ansatz of type A, which is manifestly transla- 
tional invariant in the so called uniform gauge, the gauge 
transformations of the PSG for symmetric Bosonic spin 
liquid state is found to be (we have adopted the conven- 



tion of Ref. \M) 

^Py — 'lyP x 'lxP x K 

G P = j***, 

in which T) xPx , r\ y p x = ±1, (j) x ,(p xy = 0, ir/2. 

If we require Aij to be non-zero between nearest- 
neighbor sites, the PSG should satisfy 

rixP x = -riyp x 
VxP* = ~e 21 ^. 

In the uniform gauge, Aij is only a function of Rj—Ri, 
so we can write A+ t j as Ar^^ ), in which the distance 
(d x ,d y ) = (j x - i x ,j y - i y ). By applying P x and P y 
successively, we have 

A (-d„-d y ) = (rixP m V v P m ) dm+dv &(.d m ,d y ) 
However, from the fact that A^j = — Aj.;, we have 

We thus conclude that Aj j is non-zero only between sites 
in the opposite sub-lattices. 

To show further that A k= ( Wj0 ) = 0, we need to go to 
the sub-lattice uniform gauged For (f) xy — 0, the gauge 
transformation from the uniform gauge to the sub-lattice 
uniform gauge is given by 

W. = (_!)[*=£*], 

while for <j) xy = 7r/2, it is given by 

Wi = (-l)l^, 

in which [r] means the largest integer that is not greater 
than r. In the sub-lattice uniform gauge, the pairing term 
A(d x ,d H ) has s-wave symmetry from any site in the A or 
B sub-lattice (but has opposite signs for Ajj starting 
from the A and B sub-lattices) . Thus the total contribu- 
tion to the Fourier transform of A^^ci ) from distance 
(d x ,d y ) and all the other symmetry related distances is 
proportional to 

\d x ,d y ) (cos(k x d x ) cos(kyd y ) + cos(k x d v ) cos(k y d x )) . 

Since A^ .d y ) is non-zero only when d x + d y is an odd 
integer, it is easy to see that Ai^^q) = 0. 
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